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Abstract. In a recent paper published in Nonlinear Analysis: Theory, Meth- 
ods & Applications, C. Argaez and M. Melgaard studied excited states for 
pseudo-relativistic multi-configuration methods. Their paper follows a previ- 
ous work of mine in the non-relativistic case (Arch. Rat. Mech. Anal 171, 
2004). The main results of the paper of C. Argaez and M. Melgaard are correct, 
but the proofs are both wrong and incomplete. 
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In a recent paper Q] published in Nonlinear Analysis: Theory, Methods & Ap- 
plications, C. Argaez and M. Melgaard studied the existence of excited states in 
a nonlinear model of quantum chemistry called multi- configuration. Their paper 
extends to the pseudo-relativistic setting some results which I have obtained in the 
non-relativistic case in [15] . The main results of the paper of C. Argaez and M. 
Melgaard are correct, but the proofs are both wrong and incomplete. The purpose 
of this letter is to explain why. 

In the paper of Argaez-Melgaard, there is an important confusion between Hartree- 
Fock (HF) and multi- configuration (MC) theories. Indeed, sections 5, 7 and 8 in [1] 
have simply been copied and pasted from a previous paper [7] of Enstedt-Melgaard 
on pseudo-relativistic Hartree-Fock equations, without being of any help in the 
multi-configuration case. 

Let me quickly re-explain the difference between the Hartree-Fock and multi- 
configuration models. The aim of these two nonlinear theories is to approximate 
the solutions of the many-body Schrodinger equation, describing electrons in atoms 
and molecules. The Hamiltonian of the system is the operator 



JV 

(1) ff NB l; : £(t™ l; • V:,,; 



V Xj / \xi — xA 

j=l l^k<£^N 1 1 

acting on the subspace /\^L 2 (K 3 ) of antisymmetric functions xjy) in 

L 2 ((M. 3 ) N ). The function V is the electrostatic potential induced by the nuclei 
in the system which, in the Born-Oppcnheimer approximation, are treated as clas- 
sical pointlike particles: 

M 
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V(x) :=-J2 



Here Zk > and Rk € M 3 are, respectively, the charges and positions of the nuclei. 
The total nuclear charge is Z = Yltf=i Zfe - ^ ne °P era tor T NR / R describes the kinetic 
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energy of the electrons. In non-relativistic quantum mechanics 

(2) T NR = -|. 

Relativistic particles should be described by the Dirac operator but there is no 
consistent theory for interacting systems at present. It is therefore useful to test 
some ideas on the so-called pseudo-relativistic operator 



(3) T R = oT 2 ( yl — o^A — 1 



Here a > is the (bare) fine structure constant whose inverse is the speed of light. 
Pseudo-relativistic many-body systems based on T R have been considered before in 
several important works including 15, 18 20,23,24]. The quadratic form associated 
with H K is not always bounded from below (as opposed to the case of £f NR ) but it 
is so when aZ ^ 2/ir, by the Hardy-Kato inequality. In this case H R is well defined 
by Friedrichs' method and it has similar properties as in the non-relativistic case. 
In particular, when N < Z + 1 and aZ < 2/tt, there are infinitely many eigenvalues 
Afc's below its essential spectrum [23" l l2"4 ] . corresponding to Schrodinger's equation 



(4) ff NR/R * fc = A fe * fc . 

For any normalized wave- function ty, it is convenient to define the one-particle 
density matrix 7^ by its integral kernel 



7*(x,y) = N / dx 2 --- dx N $(x, x 2 , —,x N ) i Sf(y,X2, —,x N ). 

This is a self-adjoint operator on L 2 (M. 3 ) such that ^ 7* ^ 1 and Tr(7$) = N. 
Any one-body observable can be expressed in terms of 7* only. For instance, 



V(x,)) U ) = Tr L2(R3) (T NR / R + V)^ 





where both sides are interpreted in the sense of quadratic forms. On the other 
hand, the electronic Coulomb repulsion (the second sum in formula (jTJ)) is a two- 
body operator which cannot be expressed in terms of 7$ only, except for very 
specific states like Hartree-Fock states. 

The Hartree-Fock method [T7l[2T] consists in restricting the many-body energy 
to wave-functions which are a single Slater determinant 

*hf(^1: —,xjv) = (fi A • • ■ A (p N )(x!, ...,x N ) := det(<fi(xj)), 

where (cpi,ipj) = 8{j. Such states are completely described by their one-particle 
density matrix, which is the orthogonal projection on span(y?i, ...,ipn)' 

N 
3=1 

In particular, the two-body energy (hence also the total energy) can be expressed 
in terms of 7* HF only. This fact is used in modern theoretical studies of HF-type 
models, as pioneered by Bach, Lieb, and Solovcj 2 Hi 22 . as well as in numerical 
optimization techniques [I]. 
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The multi-configuration methods are based on the observation that Slater deter- 
minants span the whole many-body space. The many-body energy is then restricted 
to wave-functions which are a finite linear combination of Slater determinants: 

(5) *Mc(afi t - )Ejv) = X! a h-i N <Ph A • • • A ifi N . 

l<ii<— <ijv<K 

The unknowns are the mixing coefficients ai 1 ,.,i N which must satisfy the normal- 
ization constraint X^l^u—wl = 1 an d the orbitals (pi,...,cpK which must be or- 
thonormal: (tfi,(pj) — 5ij- When K — N or when all the cij 1 ...j N 's vanish except 
one, we are back to the HF method. The MC energy is nothing but the many-body 
energy of such special states, expressed in terms of the mixing coefficients ai 1 ...i JV 's 
and the orbitals tpi, tptc- In general, the interaction energy cannot be expressed 
only in terms of the one-particle density matrix. The MC equations consist of a 
system of K coupled nonlinear PDE's for the tpj's, together with a (^) -dimensional 
eigenvalue equation for the a^ . ^'s. 

The existence of minimizers for non-relativistic MC was proved in a fundamental 
paper of Friesecke [9] and, later, in a paper of mine |13j . with a different method 
based on previous works by Lions [2T] and Fang-Ghoussoub [SHU] - This technique 
allowed me to construct specific critical points of the MC energy, interpreted as 
approximate excited states in a certain sense. This is what was extended to the 
pseudo-relativistic case in the paper of Argaez-Melgaard. 

Except for the sections copied from [7], the paper of Argaez-Melgaard follows 
very closely my paper [13] . I think it is fine to copy the literature, as long as the 
source is clearly mentioned, which is debatable here. The method of [13] works also 
in the pseudo-relativistic case, except for some minor steps. One difficulty is that 
the potential term involving V is not continuous for the if 1/,2 (R 3 ) weak topology: 

(6) <p n - ip weakly in H^(R 3 ) ±* I dx -» / dx. 

Jr3 f| JR3 \x\ 

Said differently, the operator (1 — A) -1 / 4 |a;| -1 (l — A) -1 / 4 is not compact. The 
result corresponding to ^ is true in the non-relativistic case in which J ff 1 / 2 (R 3 ) 
is replaced by i? 1 (R 3 ), and this was used in [13] as well as in most papers dealing 
with non-relativistic atomic models. What is true in the pseudo-relativistic case, 
however, is that the total one-body energy is weakly lower semi-continuous, that is 

(7) <p n ip weakly in H 1/2 (R 3 ) 

=>■ liminf (<p n , (T R + V)<p n ) > (ip, (T R + V)<p) 

as soon as aZ < 2/tt. This well-known fact was already employed in pseudo- 
relativistic Hartree-Fock theory by Dall'Acqua, Cstergaard S0rensen and Stock- 
meyer in [B], The reason why (J7J is true is that T R + V can be written as 
T R + V = (T R + V) + - (T R + V)_ where (T R + V) ^ and (T R + V) _ is 
compact Q It is indeed a general fact that the quadratic form of a bounded-below 
operator A is weakly lower semi-continuous for the associated topology, if and only 
if the essential spectrum of A does not go below 0. One can verify that ((TJ is 
sufficient to adapt the proof of [13] to the pseudo-relativistic case. 

4n \E\ it was even shown, using ideas of [3], that (T H + V)_ is Hilbert-Schmidt, but this is 
not necessary here. 
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In the paper of Argaez-Melgaard, the property (0 is somewhat proved in the 
appendix. It is formulated in terms of density matrices, which is fine but not 
necessary. Even though the same result in the Enstedt-Melgaard paper [7J p. 14] 
simply referred to [6], it is detailed again in the Argaez-Melgaard paper, without 
even mentioning [6]. 

What is however completely wrong in the Argaez-Melgaard paper is the treatment 
of the two-body repulsion between the electrons. This term is non- negative and the 
proof of [13] applies without any difficulty. Instead, Argaez and Melgaard "switch 
to the density operator formulation" on page 396, and they associate to any MC 
wave-function W of the form ([5]) the 'density' operator 



Then, on pages 396-397 and in the appendix, they use that the MC total energy is 
equal to the HF energy of T>. This is obviously wrong (the energy is quadratic in 
the ai t ...i N but the one-body term in the HF energy is linear). The operator T> is 
not even related to the true density matrix 7^ (except when all the a^...^ vanish 
but for one, which equals 1). 

Even if we forget for a moment that the energy cannot be expressed in terms 
of T>, most of the mathematical arguments of Argaez-Melgaard based on V are 
incorrect. The main difficulty that T> is not a non-negative operator (the a il _ iN 
can be negative), is avoided in a very questionable fashion. In several places, 
the authors seem to forget that the a,i 1 ...i N , $ have no sign and, when necessary, 
they suddenly replace ai 1 ...i N by the absolute value la^...^! (see, e.g., the liminf 
argument at the bottom of page 402). Even the nonlinear HF interaction term is 
not weakly lower semi-continuous when the density matrix has no sign, and (A. 14) 
is not correct. 

Let me end this Letter by mentioning another issue with the paper of Argaez and 
Melgaard. Their main result (Proposition 9.1) contains the statement that the limit 
of a Palais-Smale sequences for the MC energy always has a rank equal to K — 1 
or K. The non-relativistic equivalent to this statement was proved independently 
in [10l[14], based on a previous method of Le Bris [12]. It relies on the fact that 
the orbitals ipi 's of a solution to the MC equations are real-analytic away from the 
nuclear positions. This fact is not proved by Argaez and Melgaard, who instead refer 
to 'a future work' on page 397. So the corresponding statement 5 in Proposition 
9.1 is in fact not proved in the paper. Fortunately, the real-analyticity was recently 
provecH by DalPAcqua, Fournais, Cstergaard S0rensen and Stockmeyer in [5], and 
Proposition 9.1 is finally correct. 

As a conclusion, the paper of Argaez and Melgaard refers improperly to the 
literature and some parts of the proofs are wrong. 
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The proof of this in the Hartree-Fock case was made available on arXiv in the first version 
of [5], before the submission of the paper of Argaez-Melgaard. The observation that the proof 
works the same in the MC case was added in the second version of [5], after the publication of [l]. 
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